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From the Euler-Heisenberg formula we calculate the exact real part of the one-loop effective 
Lagrangian of Quantum Electrodynamics in a constant electromagnetic field, and determine its 
strong-field limit. 
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1. Introduction 

The possibility of electron-positron pair production 
from the vacuum of Quantum Electrodynamics (QED) 
vacuum was first pointed out by Sauter, Euler, Hciscn- 
berg and Schwinger Q, |^ |1| who studied the behavior of 
the Dirac vacuum in a strong external electric field. If 
the field is stronger than the critical value Ec = rn^c^ /eh, 
the energy of the vacuum can be lowered by sponta- 
neously creating an electron-positron pair. This is the 
Sauter-Euler-Heisenberg-Schwinger process. For many 
decades, both theorists and experimentalists have been 
interested in the aspects of electron-positron pair produc- 
tion from the QED vacuum by an external electromag- 
netic field. There are many reasons for the interest in 
the phenomenon of pair-production in a strong electric 
field. In addition to understand the behavior of QED in 
the strong coupling regime, the most compelling one is 
that both laboratory conditions and astrophysical events 
may now exist for observing this process. 

The above process is governed by the imaginary part of 
the Euler-Heisenberg Lagrangian, and it is this imaginary 
part which has received most attention until now. In this 
article, we supplement the existing results by exact ex- 
pressions for the real part. This determines the lowering 
of the energy of the vacuum by the external fields. 



2. QED in external electromagnetic fields 

The QED Lagrangian describing the interacting sys- 
tem of photons, electrons and positrons reads 



(1) 



where the free Lagrangians Cq^'' and Cq for electrons 
and photons are expressed in terms of Dirac field ipi^) 
and an electromagnetic field A^(a;) as follows: 



n 
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(2) 



C2 = -I- gauge-fixing term.(3) 

Here are the 4x4 Dirac matrices, ip{x) = ip^{x)'^'^, 
and Ffj,i, = d^iA^~d„A^ denotes the electromagnetic field 



tensor. Minimal coupling gives rise to the interaction 
Lagrangian 



Ant = -eiP{x)-f''iP{x)Af,{x). 



(4) 



We are using natural units h = c = 1. 

An external electromagnetic field is incorporated by 
adding to the quantum field in Q an unquantized 
external vector potential A°^, so that the total interaction 
becomes 

Ant + Ant = -e^ix)rHx) [A^ix) + Alix)] . (5) 

Quantum field theory is defined by a functional integral 
formulation for the quantum mechanical partition func- 
tion 



Z[A'']= [V%I^V'4>VA^]eyiY> 



d^x{C + Cl,,) 



(6) 



to be integrated over all fluctuating electromagnetic and 
Grassmannian electron fields. The normalized quantity 
Z[A°] gives the amplitude for the vacuum to vacuum 
transition in the presence of the external classical elec- 
tromagnetic field: 



(out,0|0,in) = 



(7) 



where |0,in) is the initial vacuum state at the time t = 
t_ —00, and (out, 0| the final vacuum state at the 
time t = t+ ^ +00. By selecting only the one-particle 
irreducible Feynman diagrams in the the perturbation 
expansion of Z[A'^] one obtains the effective action as a 
functional of A^ 



AAffK] = -iln(out,0|0,in) 



(8) 



In general, there exists no local effective Lagrangian den- 
sity A£eff whose spacetime integral is A^eS [A°] . An infi- 
nite set of derivatives would be needed, i.e., A£off would 
have the arguments A'^{x),df_iA''{x), df^d„A'^{x), . . . , con- 
taining gradients of arbitrarily high order. With 
presently available methods it is possible to calculate a 
few terms in such a gradient expansion, or a semi classi- 
cal approximation a la WKB or an arbitrary but smooth 
spacetime dependence. Under the assumption that the 
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external field A^{x) varies smoothly over a finite region 
spacetime, we may define an approximately local effective 
Lagrangian A£es[-A'^ {x)] , 

AA,sm~ I d^xAC,B[A^x)]^VAtA£,ff[A^], (9) 

where V is the spatial volume and time interval At = 
t+ 

For a large time interval At = t+ — t- —>■ oo, the 
amplitude of the vacuum to vacuum transition Q has 
the form, 



(out, 0|0, in) = e 



-i(A£o-ir/2)At 



(10) 



where AEq — £o{A°) — £o{0) is the difference between 
the vacuum energies in the presence and the absence of 
the external field, F is the vacuum decay rate, and At 
the time over which field is non zero. The probability 
that the vacuum remains as it is in the presence of the 
external classical electromagnetic field is 



|(out,0|0,in)|^ = e 



2 ^ -2ImA^off[A'' 



(11) 



This determines the decay rate of the vacuum in an ex- 
ternal electromagnetic field: 



F 2ImAAff[v4°] 



V 



AtV 



2ImA£off[A'= 



(12) 



The finite lifetime is caused by the production of electron 
and positron pairs. The external field changes the energy 
density by 



A£o ReAAes[A°] 



V 



AtV 



ReAC^siA" 



(13) 



3. Schwinger formula for pair production 



The Dirac field appears quandratically in the partition 
functional (jSJ and can be integrated out, leading to 

= J VA^ Det{i^ - e\^{x)+J^{x)] - + i^}; 

^ = rd^, (14) 

where Det denotes the functional determinant of the 
Dirac operator. Ignoring the fluctuations of the electro- 
magnetic field, the result is a functional of the external 
vector potential A'^{x): 

Z[A°] « const. X Bet{i^ ~ eJf{x) - rUe + irj}. (15) 

which is the one-loop approximation of 0{h). The in- 
finitesimal constant irj with r/ > specifies the treatment 
of singularities in energy integrals. ^From Eqs. (f7|)- (|15|l . 
the effective action Hll(l is given by 

AActfiA''] = -iTrln - ejf{x) - me + i'n] 

■ I' (16) 

ip — rrie + ir] J 



where Tr denotes the functional and Dirac trace. In phys- 
ical unit, this is of order h. The result may be expressed 
as a one-loop Feynman diagram, so that one speaks of 
a one-loop approximation. More convenient will be the 
equivalent expression 

AAff[A1 = -'-Ti-\n(^{[t^~ejr{x)]^-ml+iij} 



—d^ — m'^ + ii]J ' 



(17) 



where 



[2^ - e^(x)]2 = [^^, - eAlix)]' + ^o^-'F;,. (18) 



Using the identity 



In- 



(19) 



a2 _ ds r 
ai Jo 

Eq. H17|l becomes the Schwinger form 

/>00 J 

AA^aW] = 7; ^^~^s(ral-^n) (2O) 

2 Jo ^ 

■ Tr(a;|e^H[*a.--^;i(-)l'+i-"'^f^;:4 - e'^^^'l^). 



where a^^ = 



d^^Al 



d,yA°^ , and 



{x\{- ■ -jla;) are the diagonal matrix elements in the local 
basis \x). This is defined d la Dirac by the matrix el- 
ements with the momentum eigenstates |A;) being plane 
waves: {x\k) — e^**^^. The symbol Tr denotes integral 
J d'^x in spacetime and the trace in spinor space. For con- 
stant electromagnetic fields, the integrand in (|20|l does 
not depend on x, and cr^^i^^^ commutes with all other 
operators. This will allow us to calculate the exponential 
in Eq. 1)20(1 explicitly. The presence of —irj in the mass 
term ensures the convergence of integral for s — > oo. 

If only a constant electric field E is present, it may be 
assumed to point along the z-axis, and one can choose a 
gauge such that = —Et is the only nonzero compo- 



nent of A° . Then one finds 



trexpis 



-a'^'F" 

2 ''r 



~ 4cosh(sei?), 



(21) 



where the symbol tr denotes the trace in spinor space. 
Using commutation relation [c^oiS;*^] = 1, where x° = t, 
one computes the exponential term in the effective action 
(EOl (c.e.g. 0) 



{x\ exp IS 
eE 



{id^-eAl{x))' + -a^-^F^^, \x) 



(27r)2is 



coth(ei?s). 



(22) 



The second term in Eq. (|20|) is obtained by setting E ~ 
in Eq. (|22|l . so that then the effective action (|20|l yields. 



1 r pGo J 
AAff = j d\ ^ [eEscot\i{eEs) - 1] 



-2s(mg —ir}) 



(23) 
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Since the field is constant, the integral over x gives a 
volume factor, and so the effective action (|2U|) can be 
attributed to the spacetime integral over an effective La- 
grangian © 



1 /"Oo ^ 

(24) 

By expanding the integrand in Eq. 1241) in powers of e, 
one obtains, 



4r [eEs coth(e£:s) - 1] e" 



s(r] 



~is{m^—iri) 



_3s 45 
The small-s divergence in the integrand, 
1 



3 2(27r)2 



-is(m^-ir]) 



(25) 



(26) 



is proportional to the electric field in the original Maxwell 
Lagrangian. The divergent term H2t)|l can therefore be 
removed by a renormalization of the field E. Thus, we 
add subtract counterterm in Eq. (|24|l and form, 



AC 



1 



off 



2(27r) 



ds 



eEscoth{eEs) - 1 - y-E'^s^ 



-is(m^—i7]) 



(27) 



Remembering Eq. (|12|l . we find from (|27|l the decay rate 
of the vacuum per unit volume 



r 

V 



(2^)2 



-Im 







ds 
Is 



eEscoth{eEs) - 1 - y^^^s^ 



-is{m^ —11]) 



(28) 



The integral H28|) can be evaluated analytically by pro- 
ceeds by the method of residues. Since the integrand is 
even, the integral can be extended to the entire s-axis. 
After this, the integration contour is deformed to enclose 
the negative imaginary axis and to pick up the contri- 
butions of the poles of the coth function at s = mr/eE. 
The result is 



r 

V 



aE^ ^ 1 



-^exp 



niiEc 
E 



(29) 



where Ec — rr? je is the value of critical field. This result 
due to Schwinger Q is valid to lowest order in h for 
arbitrary constant electric field strength. 

An analogous calculation for a charged scalar field 
yields 



r 

V 



aE' 
2^ 



2.— ;^exp 



rniEr 



(30) 



These Schwinger results complete the development. The 
leading n = 1 -terms agree with the WKB results ob- 
tained by Sauter and Heisenberg-Euler Q. 



4- Pair production in constant electromagnetic fields 

If the constant external field has both E and B 
nonzero, the effective action will be a function of the 
two Lorentz invariants. We can now go to an arbitrary 
Lorentz frame by expressing the result in terms of the 
scalar S and the pseudoscalar P, that can be formed 
from arbitrary electromagnetic fields B and E fields: 



1 



1 



P 



If PA"' 



E B, 



(31) 

where F'^'^ = e'^^^'^Fx^. is the dual field tensor. It is useful 
to define the related invariants e and (3 as the solutions 
of the invariant equations 



£2 - /32 = e2 - B^ = 25-, e/3 = E • B = P, 
which read explicitly 



(32) 



1 



= -^y v/(E2 - B2)2 + 4(E • B)2 ± (E2 - B2) 

V ^ 



(52+P2)l/2±^, 



(33) 



Then exists a special Lorentz frame to be called the 
center- of- fields frame, where the fields Bcf and Ecp are 
parallel. In this frame, /3 — |Bcf| and e — |Ecf|- 

Relations (|32|l are invariant under the discrete duality 
transformation: 



IBI 



-i\E\ 



lEl 



«|B| 



i.e., under 



f3 



-le, e ■ 



if3. 



(34) 



(35) 



This implies that effective action for a pure magnetic 
field can be simply obtained by replacing E\ — > iB in 
Eqs. (|22I27|I . so that 



and 



{x\ exp IS 

- e-^ 
~ {2nyis 



1 



eAl{x)f + -a^'^F;^ 



cot{eBs), 



(36) 



2(27r)2 

— is(m^ 



ds 



eBscot{eBs) — 1 



-irj) 



(37) 



If both electric and magnetic fields E and B are 
present, we assume the parallel Ecf and Bcf fields point 
along the z-axis in center- of -fields frame. We further 
choose a gauge such that only A° = —Ecpt, Ay — Bcfx^ 
are nonzero. The exponential in the effective action 
Eq. H20() can then be factorized into a product of the 
magnetic part and the electric part. Based on the dual- 
ity, we obtain directly 



tr exp is 



-af"'F'' 
2 



4cosh(sei?cF) cos(sei3cF), (38) 



4 



and 



Effective nonlinear Lagrangian 



{x\ expis {[id^ 
1 eEcF 



{2ny 



IS 



coth{seEcF)^^^cot{seBcF)- (39) 



Performing the same substraction as before, we obtain 
the effective Lagrangian 



AC. 



1 



cff 



2(27r)2 



ds r 



e^-EcF-Bcps^ coth(sei5cF) 



cot(seBcF) - 1 - y(-^CF 



BiF>' 



s(ri 



irj) 



(40) 



In an arbitrary Lorentz frame, Eqf and Bqf are replaced 
by the invariants e and 



AC. 



cff 



2(27r)2 

1 



ds r 
3 



e^e/3s^ coth(ees) cot(e/3s) 



(41) 



The decay rate becomes 
1 



r 

V 



-Im 



ds ' 

73 



e^e/3s^ coth(ees) cot{e(3s) 



l-y(£'-/3V 



(42) 



The integral is evaluated as in Eq. (|29f) by the method 
of residues, and yields Q 



r 

V 



OLE X ^ 

r2 



TT^ tanhn7r/3/e 

71—1 ' 



exp 



nnEr 



(43) 



which reduces for /3 (B = 0) correctly to (|29|l . 
The analogous result for bosonic fields is 



ae 



r 



V 27r2 ^ 

n— 1 



(-1)" n7r/3/£ 

sinh nnP/s 



exp 



rnrEc 



(44) 



Note that the magnetic field produces in the fermionic 
case a extra factor {mr P / e) / tanh.{mT f3 / e) > 1 in each 
term which enhances the decay rate. The bosonic series 
(|44|l . on the other hand, carries in each term a suppres- 
sion factor {mrP/e)/ sinhmr/S/e < 1. 

The decay rate T/V gives the number of electron- 
positron pairs produced per unit volume. The prefac- 
tor can be estimated on dimensional grounds. It has 
the dimension of E'^/fi, i.e., m^c'/h'^. This arises from 
the energy of a pair 2meC^ divided by the volume whose 
diameter is the Compton wavelength h/rrieC produced 
within a Compton time h/nieC^- The exponential factor 
suppresses pair production as long as the electric field 
is much smaller than the critical electric field Ec- The 
general results H43|) , H44|l was first obtained by Schwinger 
[Sj for scalar and spinor electrodynamics. 



In this section, we evaluate further the effective La- 
grangian (|41|l for arbitrary constant electromagnetic 
fields E and B. Making the expansions, 



ees coth(ees) = 



e/9s cot(e/3s) 



oo 

E 

n— — oo 
oc 

E ■ 



(s^ + r^y 

o2 



m— — oo 

we obtain the finite effective Lagrangian, 



AC 



1 



off 



2(27r) 



E ' 



n,m— — oo 



ds- 



mr/ee, (45) 
E m7r/e/3,(46) 



-is{7n'^ —ir]) 



{S^+T^)i 

where divergent terms ny^0,TO = 0,n = 0,r7i^0 and 
n = TO = are excluded from the sum, as indicated by 
a prime. The symbol 5ij = 1 — 5ij denotes the com- 
plimentary Kronecker-(5 which vanishes for i = j. The 
divergent term with n = to = is eliminated by the zero- 
field subtraction in Eq. 141(1 . while the divergent terms 
n 0,m — and n = 0, to 7^ 




(48) 



are eliminated by the second subtraction in Eq. 141|l . 
This can be seen by performing the sums 

k 



h (?) «'* s 



1 



e/5 



C(fc), (49) 



where C(fc) = l/n*' is the Riemann function. 

The infinitesimal irj accompanying the mass term in 
the s-integral (|47l) is equivalent to replacing e-»s(m;-«»?) 
by e-»s(i-«'?)™e . This implies that s is to be inte- 
grated slightly below (above) the real axis for s > 
(s < 0). Equivalently one may shift the (— t,„) vari- 
ables slightly upwards (downwards) to Tm + iiJ {—Tm — ir]) 
in the complex plane. In order to calculate the finite ef- 
fective Lagrangian (|47|l . the factor e-«s(i-«'/)™e is divided 
into its sin and cos parts: 



4(27r)2 



E ' 



n.m— — QO 



sds 



->mO 



sin[s(l — iri)m^ 



eJ I 



(50) 



A /-cos 

^^off 



00 



2(2^)2 — ^0 



sds 



JmO 



(s2 -t- r2 



cos[s(l — irj)ni\ 



(51) 
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The sin part H50|l has an even integrand allowing for an 
extension of the s-integral over the entire s-axis. The con- 
tours of integration can then be closed by infinite semicir- 
cles in the half-plane, the integration receives contribu- 
tions from poles ±t„i, ±iT„, so that the residue theorem 
leads to 



oft 



— \ - coth I — - 1 exp{-n7rEc/e)i52) 



—I 



^2 $Z 

m— 1 



2(27r)- m 



coth 



e 

/ rmrs 



exp{irmrEc/(3){53) 



The first part H52|l produces the exact non-perturbative 
Schwinger rate l|43|) for pair production. 

Shifting s ^ s~ irj, we rewrite the cos part of effective 
Lagrangian H51|l as 



2(27r)2 



n,m — — ck; 



, cosjsml) 

'n ^ 'm 



sSnO 



irj s-^ + T„ — irj 



(54) 



In the first term of magnetic part, singularities s — 
Tm, (to > 0) and s = — r^, (m < 0) appear in integrating 
s-axis. We decomposes, 



i-d(s - Tm) + *-^'^(s + '^m) 



(55) 



where 7^ indicates the principle value under the integral. 
The integrals over the (5-functions give 

^5^cff = *7^S y. 3 coth f exp(*m^i?,//3). 



2(27r)2 ^ TO 

^ ' m=l 



(56) 

which exactly cancel the second part (|53|l of the sin part 

A /-sin 

It remains to find the principal-value integrals in 
Eq. ini, 



(A/:s^)p 



2(2^)- — T. 



2 

n 

Sf^mO S(5„o 



1 /"^ 

^"P / dscos(smg) 

(57) 



We rewrite the cos function as cos(sTOg) — (e**™= + 
g-2sm^-j^2 and make the rotations of integration contours 
by ±7r/2 respectively, 

1 1 /-oo 



Snare 
■ - (T„m2) 



(58) 



Using the formulas (see 3.354, 8.211.1 and 8.211.2 in Ref. 
|6jj) 



J{z) = V ds- 
Jq 



se 



e~^Ei(z) + e'^Ei(-z) , 
(59) 



where Ei(z) is the exponential-integral function. 



Ei(z) 



' k=l 



(60) 



we obtain the principal- value integrals H58|) . 

1 ^ , 1 



2(27r)2 



E ' 

n.m— — oo 
,2 



' m ' ' r) 



SmoJ{iT,nmJ - SnoJiTninJ . (61) 



With the help of the series and asymptotic represen- 
tation (see formula 8.215 in Ref. 6J) of the exponential- 
integral function Ei(z) for large z, corresponding to weak 
electromagnetic fields (e ^ 1, f3 <^ 1), 



J(z) 



1 6 

T2 ~ 74 



z' 

and Eq. we find. 



120 5040 

~~6 



362880 



r8 



(62) 



2(27r) 

'SrnO 



oo 



n,m— — OO 

1 6 



•2to4 
1 



120 
120 



m e 'me 



}(63) 



Applying the summation formulas H49|l . the weak- field 
expansion H63|l is seen to agree with the Heisenberg and 
Euler effective Lagrangian 



2a' 
A5mt 
64™^ 



315m; 



{(E2_b2)2 + 7(E-B)2} 
{2(e2-B2)3 + 13(e2-B2)(E • Bf} 



+ (64) 

which is expressed in terms of a powers series of week 
electromagnetic fields up to 0{a^). The expansion co- 
efficients of the terms of order n have the general form 
m'^/{Ec)"'. As long as the fields are much smaller than 
Ec, the expansion converges. 

On the other hand, using the series and asymptotic 
representation (formulas 8.214.1 and 8.214.2 in Ref. 0) 
of the exponential-integral function Ei(z) for small z ^ 
1, corresponding to strong electromagnetic fields (e ^ 
1, /?»!), 



J(z) 



eMn(z) + e-Mn(-z) +7 + 0(2), (65) 
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with 7 = 0.577216 being the Euler-Mascheroni constant, 
we obtain the leading terms in the strong-field expansion 
of Eq. (jnU, 



E 



^2 ^ ^2 

n,m— — oo ""^^ 



2(27r)2 

l„o ln(T„mg)-5™o ln(T„mg) 



.(66) 



In the case of vanishing magnetic field B = and m = Q 
in Eq. (|66|l . we have 



1 °° 1 



n— 1 



(67) 



for a strong electric field E. In the case of vanishing 
electric field E = and n = in Eq. we obtain for 
a strong magnetic field B, 



1 °° 1 

^ ' m=l 



!^ V — 



In 



B 



,(68) 



a result first considered by Weisskopf 



6. Concluding remarks 



We have presented in Eqs. H52I53I56I61|) closed form 
result for the one-loop eS'ective Lagrangian A£cff 
for arbitrary strength of constant electromagnetic fields. 
The results will receive fluctuation corrections from 
higher loop diagrams. These carry one or more factors a, 
a^, ... and are thus suppressed by factors 1/137. Thus 
results are valid for all field strenghts with an error no 
larger than roughly 1%. If we include, for example, the 



two-loop correction, the first term in the Heisenberg and 
Euler effective Lagrangian (|64|l becomes 



2a^ 



Admji ( 
7 ( 1 ''''^ 



1 + ^) (E^-B^^^ 
97r 



2527r 



(E-B)2|. 



(69) 



Readers can consult the recent review article 9j, where 
one finds discussions and computations of the effective 
Lagrangian at tow-loop levels, in homogeneous and inho- 
mogeneous fields. 

The interaction of an external electromagnetic field 
with the QED vacuum leads to the appearance of nonlin- 
ear real and imaginary parts in the effective Lagrangian 
A£off of the electromagnetic field. Electron and positron 
pair production is only one aspect of the phenomena con- 
necting to these nonlinearities. As we remarked before 
the rate T/V of pair production (|29|l is exponentially 
small when E ^ Ec- 

Unfortunately, it seems inconceivable to produce a 
macroscopic static field with electric field strengths of the 
order of the critical value in any ground laboratory to 
directly observe the Sauter-Euler-Heisenberg-Schwinger 
process of electron-positron pair production in vacuum. 
In Ref . llOl ■ we discuss some ideas: (i) heavy-ion col- 
lisions |ll|. (ii) the focus of coherent laser beams [l^l 
and (iii) electron beam-laser collisions |l3j | , to experimen- 
tally create a transient electric field E E(, 'va earth- 
bound laboratories, whose lifetime is expected to be long 
enough (larger than Ti/rrieC^) for the pair production pro- 
cess to take place. In Ref. 0, we also discuss (i) the 
idea that the supercritical conditions can be reached on 
the time scale of gravitational collapse (GM/c^, 10~^sec. 
for M = IOM0) around an astrophysical object 0| and 
(ii) the application of Sauter-Euler-Heisenberg-Schwinger 
process to the black hole physics in connection with the 
observed astrophysical phenomenon of gamma ray bursts 

m 
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